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ABSTRACT 

We examine the flaring behaviour of the cataclysmic variable AE Aqr in the 
context of the ‘magnetic propeller’ model for this system. The flares are thought 
to arise from collisions between high density regions in the material expelled from 
the system after interaction with the rapidly rotating magnetosphere of the white 
dwarf. We calculate the first quantitative models for the flaring and calculate the 
time-dependent emergent optical spectra from the resulting hot, expanding ball of 
gas. We compare the results under different assumptions to observations and derive 
values for the mass, lengthscale and temperature of the material involved in the flare. 
We see that the fits suggest that the secondary star in this system has Population II 
composition. 

Key words: accretion, radiative transfer, stars: flares, stars: individual: AE Aqr, 
novae, cataclysmic variables 


1 INTRODUCTION 

AE Aqr is an unusual cataclysmic variable star exhibiting 
bizarre phenomena that can now be interpreted in the con¬ 
text of a mag netic propeller that throws g as out of the bi¬ 
nary system (Wynn, King & Horne 1997). The magnetic 
propeller efficiently extracts energy and angular momentum 
from the white dwarf and transports it via magnetic fields 
to gas stream material which it ejects from the binary svs- 


tem. D etailed understanding of the magnetic propeller ef- 
fects may therefore help us to unlock some of the secrets of 
magnetic viscosity in accretion flows. 

In the AE Aqr binary sys tem, a slightly evolved K star 
(Welsh, Horne & Corner 1995) that overflows its Roche lobe 
is locked in a 9.88 hour orbit with a rapidly spinni ng magne¬ 
tized white dwarf. Coherent oscillations in optical ( Patterson 
197E ), ultravio let ( jSracleous et al. 1994 ), and X-ray (Patter- 
5on et qL 1981)1) lielitcurvos reveal the white dwarf’s 33s 


spin 


period. I The oscillation has two unequal peaks per spin cv- 
cle, consistent with broad hotspots on opposite si des of the 
white dwarf ~ 15° above and below the equator (Eracleous 
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et al. 1994). These oscillations are strongest in the ultravio¬ 
let, where their spectra show a blue continuum with broad 
Ly q absorption consistent with a white dwarf (log q « 8) 
atmosphere with T ~ 3 x 10 4 K (Eracleous et al. 1994). The 
simplest interpretation is accretion heating near the poles of 
a magnetic dipole field tipped almost perpendicular to the 
rotation axis. 

_ An 11-year study of the optical oscillation period jdq 


Jager 1994) shows that the white dwarf is spinning down 


at an alarming rate. Something extracts rotational energy 
from the white dwarf at a rate Ilou) ~ 60i >L V ie. some 60 
times the luminosity of the system. We now believe that to 
be a magnetic propeller. Th is m odel was first proposed by 
Wynn, King fe Horne (1995) and Eracleous fe Hor ne (1990 


and expanded upon in Wynn, King fe Horne (1997| ) with 
comparison of observed and modelled tom ograms. Further 
work on the fl arin g region was reported in Welsh, Horne & 


Comer (199g| ) and Horne (199£). 


The gas stream emerging from the companion star 
through the LI nozzle encounters a rapidly spinning mag¬ 
netosphere. The rapid spin makes the effective ram pressure 
so high that only a low-density fringe of material becomes 
threaded onto field lines. Most of the stream material re¬ 
mains diamagnetic, and is dragged toward co-rotation with 
the magnetosphere. As this occurs outside the co-rotation 
radius, this magnetic drag propels material forward, boost¬ 
ing its velocity up to and beyond escape velocity. The ma¬ 
terial emerges from the magnetosphere and sails out of the 
binary system. This process efficiently extracts energy and 
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angular momentum from the white dwarf, transferring it via 
the long-range magnetic field to the stream material, which 
is expelled from the system. The ejected outflow consists of 
a broad equatorial fan of material launched over a range of 
azimuths on the side away from the K star. 

The material stripped from the gas stream and threaded 
by the field lines has a different fate, one which we be¬ 
lieve gives rise to the radio and X-ray emission. This ma¬ 
terial co-rotates with the magnetosphere while accelerating 
along field lines either toward or away from the white dwarf 
under the influences of gravity and centrifugal forces. The 
small fraction of the total mass transfer that leaks below 
the co-rotation radius at ~ 5 R „a accretes down field lines 
producing the surface hotspots responsible for the 33s os¬ 
cillations. Particles outside co-rotation remain trapped long 
enough to accelerate up to relativistic energies through mag¬ 
netic pumping, eventually reaching a sufficient energy den- 


sity to b reak away from the magnetosphere (Kufjpersetal 
1997 |). j rhe resulting ejection of balls of relativistic magne¬ 
tized plasma is thought to give rise to the f l aring radio emis- 
sion ( iBastian^JPulk fe Chanmugam 1988a ; Bastian, Dulk & 
Chanmigam 1988b). 


This paper addresses the optical and ultraviolet vari¬ 
ability seen in AE Aqr. In many studies the lightcurves ex¬ 
hibit dramatic flares, with 1-10 minute rise and fall times 


(Patterson 1979; van Paradijs ; Kraakman fc van Ameron- 


gen 19$9|: p3ruchl99i ; Weish"^ome‘‘fc*‘Oke‘T993| ). The 
flares seem to come in clusters or avalanches of many super¬ 
imposed individual flares separated by quiet in tervals of 
gradua lly declin i ng line and cont inuum emission (Eracleous 
& Horne 1996; Patterson 1979| ). These quiet and flaring 
states typically last a few hours. Power spectra computed 
from the lightcurves have a po wer-law form, with large r 
amplitudes on longer timescales. Elsworth & James (1982) 
found that the index in A oc v a was -1. That is, that the am¬ 
plitude of the fla re or flicker v aries inversely as the frequency 
of its occurence. Bruch (1992) examined several datasets and 


found values for a in the range —0.71—1.64. Such power-law 
spectra are often associated with physical processes involv¬ 
ing self-organized critic ality, for ex ample earthquakes, snow 
or sandpile avalanches (Bak 1996). Similar red noise power 
spectra are seen in active galaxies, X-ray binaries, and other 
cataclysmic variables, and is therefore regarded as charac¬ 
teristic of accreting sources in general. However, flickering 
in other c ataclysmic v ariables typically has an amplitude 
of 5-20% (Bruch 1992), contrasting with factors of several 
in AE Aqr. If the mechanism is the same, then it must be 


weaker or dramatically diluted in other systems. 


The optical and ultraviolet spectra of the AE Aqr flares 
are not understood at present except in the most general 
terms. The lines and continua rise and fall together, with lit¬ 
tle ch ange in the equivalent w idths or ratios of the emission 
lines ( Eracleous &: Horne 1990 ). This suggests that the flares 
represent changes in the amount of material involved more 
than changes in physical conditions. The Balmer emission 
lines decay somewhat more slowly than the optical contin- 
uum - perhaps reveali ng a recombination time delay ( |Wclsh, 
Horne jz Corner 1998). Ultraviolet spectra from HST reveal 


a wide range of lines representing a diverse mix of ionization 
states and densities. Eracleous & Horne (1996) conclude: 


lines we suggest that in a large fraction of the line-emitting gas 
the density is in the range n ~ 10 9 10 11 cm -3 . It is likely that 
denser regions also exist.” 


thus setting a lower limit on the density range in the flar¬ 
ing region. Such spectra suggest shocks. The CIV emission 
is unusually weak, suggesting non-solar abundances. For ex¬ 
ample, carbon depletion may occur if CNO-processed ma¬ 
terial is being transferred from the secondary star, which is 
an evolved star being whittled down by Roche lobe over¬ 
flow. So far no quantitative fits to the spectra have been 
achieved either by shock or photo-ionization models. IUE 
observations by Jameson, King & Sherrington (1980) de¬ 
rive an emission measure (VmifVe) from Lya of 10 bi m~ 3 . 
UBVRI colour photometry of flares by Beskrovnaya et al.( 
1996) give colour indices close to those for a blackbody in 
the 15 000-20 000K range with an emitting area ~ 10 16 m 2 . 


What mechanism triggers these dramatic optical and 
ultraviolet flares? Clues come from multi-wavelength co¬ 
variability and orbital kinematics. Simultaneous VLA and 
optical observations show that the radio flux variations oc¬ 
cur on similar timescales but are not correlated with the 
optical and ultravi olet flares, which therefore require a dif¬ 
ferent mechanism (Abada-Simon et al. 1995). It was pro¬ 
posed that the flares represent modulations of the accretion 
rate onto the white dwarf, so that they should be correlated 
with X-ray variability. Some correlation was found, but the 
correlation is not high. However, HST observations discard 
this model, because the ultraviolet oscillation amplitude is 
unmoved by transitions be tween the quiet and flaring states 
(Eracleous & Horne 1996). This disconnects the origins of 


the oscillations and flares, and the oscillations arise from 
accretion onto the white dwarf, so the flares must arise else¬ 
where. 

Further clues come from emission line kinematics. The 
emission line profiles may be roughly described as broad 
Gaussians with widths ~ 1 000 km s~ J , though they of¬ 
ten exhibit kinks and so metimes multiple peaks. Deta iled 
study of the Balmer lines (Welsh, Horne & Corner 1998) in¬ 
dicates that the new light appearing during a flare can have 
emission lines shifted from the line centroid and somewhat 
narrower, ~ 300 km s _1 . Individual flares therefore occupy 
only a subset of the entire emission-line region. 

The emission line centroid velocities vary sinusoidally 
with orbital phase, with semi-amplitudes ~ 200 kin s -1 and 


maxi mum redsliift near phase ~0.8 (Welsh, Horne & Corner 


199S). These unusual orbital kinematics are shared by both 


“Based on the critical densities of the observed semiforbidden 


ultraviolet and optical emission lines. The implication is that 
the flares arise from gas moving with a ~ 200 km s -1 velocity 
vector that rotates with the binary and points away the 
observer at phase 0.8. This is hard to understand in the 
standard model of a cataclysmic variable star, though many 
cataclysmic variables show similar a nomalous emission-line 
kinematics (Thorstensen et al. 1991). Kepler velocities in an 
AE Aqr accretion disc would be > 600 km s _1 (though we 
believe no disc to be present.) The gas stream has a similar 
direction but its velocity is ~ 1 000 km s -1 . A success of 
the magnetic propeller models is its ability to account for 
the anomalous emission-line kinematics. The correct velocity 
amplitude and direction occurs in the exit fan just outside 
the Roche lobe of the white dwarf. But the question remains 
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Fireball Models for Flares in AE Aquarii 3 


of why the flares are ignited here, several hours after the gas 
slips silently through the magnetosphere. 

The key insight which solved this puzzle was the realiza¬ 
tion that the magnetic propeller acts as a blob sorter. More 
compact, denser, diamagnetic blobs are less affected by mag¬ 
netic drag. They punch deeper into the magnetosphere and 
emerge at a larger azimuth with a smaller terminal velocity. 


P = 


M 

QV 


where the stream cross-section 


Q « 2.4 x 10 1 


( Ta ) 

VlO 4 KJ 


P orb (h) m 


( 3 ) 


( 4 ) 


has been derived by several authors (Papaloizou fc Bath 


These c 

ompact blobs can therefore be overtaken by ‘fluffier’ 

1975; Meyer & Meyer-Hofmeister 1983 

; IHameury, King & 

blobs e. 

ected with a larger terminal velocity in the same di- 

Lasota 1981; Ritter 1988 

Kovetz, Prialnik & Shara 1988; 


rection, having left the companion st ar somewhat later but 
spe nt le ss time in the magnetosphere ( Wynn, King fe Horne 
1997j)-Jhe result is a collision between two gas blobs, which 
can give rise to shocks and flares. Calculations of the tra¬ 
jectories of magnetically propelled diamagnetic blobs with 
different drag coefficients indicate that they cross in an arc¬ 
shaped region of the exit stream, in just the right place 
t o account for the orbital kinematics of t he emission lines 
(Welsh, Horne & Comer 1998; Horne 1999). Figure [Ij shows 
how these trajectories map to a locus of points in the lower 
left quadrant of the Doppler map that is otherwise difficult 
to populate. 

There remains the problem of developing a quantita¬ 
tive understanding of the unusual emission-line spectra of 
the blobs. We have two aims for this paper: the principle 
aim is to see if the observed peak optical flare spectrum can 
be reproduced by conditions appropriate to the aftermath 
of a collision between blobs. The secondary objective is to 
develop as far as possible an understanding of the mecha¬ 
nisms by which the observed spectra are formed and evolve. 
We have approached both goals in the spirit of attempting 
to find the simplest models that explain the observations. 

In section ^ we outline the basic assumptions used in 
our models. We go on in section [ 3 ] to consider the radia¬ 
tive transfer problem and calculate analytic expressions for 
the behaviour of the optical flare spectra and lightcurves 
and outline some numerical considerations . In section we 
present the results of numerical simulations of the flare be¬ 
haviour, both lightcurves and spectra, and consider further 
the limits of applicability of our method. We discuss the in¬ 
terpretation of our results in section ^ in terms of the fireball 
model. Finally we summarise our results in section 


2 FIREBALL MODELS 
2.1 Rough Estimates 

We derive rough estimates for the physical parameters asso¬ 
ciated with the colliding blobs by using the typical rise time 
t ~ 300 s and an observed optical flux /„ ~ 50 mjy with a 
closing velocity V ~ 300 km s -1 for the blobs. We estimate 
a mass transfer rate 


■ 2 uL v 2 vf v d 2 14 , —i 

M ~ = ——— « 10 kg s 


(1) 


V 2 V 2 

using a distance of lOOpc. This compares to a value from 
the standard evolutionary equation 

5 

-Ml S3 6 x lO -10 (y) 3 = 3 x 10 14 kg s' 1 (2) 


( Frank, King fc Raine 1992|) . The mean density of the over¬ 
flowing gas stream in a CV can be calculated from 


T s = 4 500 K (Skidmore et al. 2002), we have a mean density, 


p ~ 3.2 x ( 300 k ' m!r , ) I*?™- 


t«7lH 


= 3.2 x 10 19 m -3 


( 5 ) 


( 6 ) 


where fx ss 0.6 is the mean molecular weight appropriate for 
a fully ionized gas with solar composition. This is consistent 
with the lower limits from HST uv observations discussed 
earlier. It must be remembered that this is the mean density 
in a smooth stream. A stream composed of discrete blobs of 
material would have range of densities about this value. 

We estimate the total mass involved in the collision to 
be 


M ~ Mt ~3x 10 16 kg. 


( 7 ) 


The typical pre-collision lengthscale a of the problem is given 

by 


a ~ Vt ~ 9 x 10' m. 


( 8 ) 


Following the collision we might expect a somewhat smaller 
lengthscale eg.~ 5 x 10 7 m. The above density then suggests 
an emission measure Qan 2 « 10 63 nW 3 consistent with the 
above IUE measurements for an optically thick Lya line. 

The initial post collision temperature T following a 
strong shock 


T ■ 


3 / urrinV 2 
16 k 


1.2 x 10° K. 


( 9 ) 


Finally we have the energy involved in the collision and max¬ 
imum possible flare energy E 


E ' 


MV 2 


3 x 10 26 J. 


( 10 ) 


A summary of typical values is provided in Table hi. 


2.2 Initial Conditions 

We envision an expanding fireball emerging from the after- 
math of a collision between two masses mi + m 2 = M. In 
the centre of mass frame, these masses approach with initial 
velocities satisfying m-ivi = m. 2 V 2 - The initial kinetic energy 
in the centre of mass frame is 


E = 


mivf + m2«2 


MV 2 q 

2(1 + qf 


( 11 ) 


where V = vi + V2 is the (frame-independent) relative ve¬ 
locity and q = mi/m,2(= V2/V1) is the mass ratio. Note that 
for a given M and V a range of energies is possible, ranging 
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P— 9.88h M, —0.9 M 2 — 0.57 K, = 1 09 K 2 =173 i = 60 




Figure 1. Trajectories of diamagnetic blobs passing though the AE Aquarii system with the corresponding Doppler tomogram. 
The open circles mark the time of flight in units of 0.1 P or b- Asterisks and grey patches mark the locations where lower density 
blobs overtake and collide with more compact blobs, producing fireballs. This occurs in the lower left quadrant at the location 
corresponding to the observed emission lines. 


Quantity 

Value 

Observed Flux f v 

50 

mjy 

Closing Velocity V 

300 

km s -1 

Flare Risetime £ r i se 

300 

s 

Mass Transfer Rate M 

10 14 

kgs- 1 

Fireball Mass M 

3 x 10 16 

kg 

Pre-Collision Lengthscale a 

9 x 10 7 

m 

Initial Temperature T 

10 6 

K 

Total Energy E 

3 x 10 26 

J 

Typical Density p 

3 x 10 -8 

kg m -3 

Number Density n 

3 x 10 19 

m -3 

Column Density N 

3 x 10 27 

m -2 


Table 1 . Estimates of typical values for flare quantities. 


from zero for very extreme mass ratios up to a maximum of 
Eq = MV 2 /8 for equal masses. 

As the collision progresses, some of the energy E is con¬ 
verted to thermal energy through dissipation in shocks. The 
initial dissipative phase lasts a short time 


ti 


(q 

V 


= 300 


(wm)( 


V 


300 km s 


( 12 ) 


The ratio E/M of energy per mass sets the initial tem¬ 
perature Ti and sound speed c s ,i through 


E _ -ykTi 
M prrtH 


(13) 


where a specific heat ratio 7 = 5/3 is appropriate for fully 
ionized atomic gas. For V ~ 300 km s -1 , Ti ~ 10 6 K, so that 
the initial ball of hot gas is atomic and completely ionized. 

Since V ~ c s ,i, sound waves have time to cross the initial 
fireball during the collision time ti, and we therefore expect 
a roughly uniform temperature profile inside a,. 


2.3 Expansion 


With nothing to hold it back, the hot ball of gas expands 
at the initial sound speed, launching a fireball. We adopt a 
uniform, spherically symmetric, Hubble-like expansion V = 
Hro = Hr/(3, in which the Eulerian radial coordinate r of 
a gas element is given in terms of its initial position r 0 and 
time t by 

r(ro,t) = r 0 + v(r 0 )t = ro + Hr 0 t = r 0 (l + Ht ) = r 0 (3. (14) 


This defines an expansion factor (3 which we can use as a 
dimensionless time parameter. The ‘Hubble’ constant is set 
by the initial conditions Hai « v{a/) « c s ,i- 

Uniform free expansion is a suitable approximation be¬ 
cause the flow becomes supersonic as it expands and cools. 
If we can determine parameters at some time t = 0 ie. (3 = 1 
for the lengthscale (a 0 ), temperature (To) and mass (M) we 
can derive the time evolution as follows. We adopt a Gaus¬ 
sian density profile 

p(r,t) = p 0 (3~ 3 e~ v2 (15) 


where 

M 



is the initial central density and 
_ ro _ r 
^ ao a 


(16) 


(17) 


is a dimensionless radius coordinate, the radius r scaled to 
the lengthscale a = (3a 0 . 

Although this Gaussian density profile is only a guess, 
it is motivated by the Gaussian shapes of observed velocity 
profiles, and by the thought that the initial thermal velocity 
distribution in the hot gas will map into the density profile 
of the expanding fireball because the fast particles travel 
farther than slow ones. 
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2.4 Cooling 

We consider three cooling schemes: adiabatic, isothermal 
and radiative. The adiabatic fireball cools purely as a re¬ 
sult of its expansion and corresponds to a situation where 
the radiative and recombination cooling rates are negligible. 
In contrast, the isothermal model maintains a fixed temper¬ 
ature throughout its evolution. A truly isothermal fireball 
would require a finely balanced energy source to counter¬ 
act the expansion cooling. However, it may an appropriate 
approximation to a situation where a photospheric region 
dominates the emission and presents a fixed effective tem¬ 
perature to the observer as a result of the stong dependence 
of opacity on temperature. The radiative models cool adia- 
batically throughout and as a result of radiation from a thin 
zone near the photosphere which we model by immediately 
dropping the temperature to T p ~ 10 4 K at this surface. 


2-4-1 Adiabatic Cooling 

During adiabatic expansion, P oc p 1 and T oc P/p oc p 7 ^ 1 . 
If we consider the evolution of a given gas element then, 
with p oc /3 -3 , the temperature decreases with time as 

T(r,t) =T 0 (ro)/3“ 3(7 “ 1) . (18) 


We see that an initially uniform temperature distribution 
remains uniform and so, for this case, we can write, more 
simply, 

T(r,t) =T 0 /r 3(7-1) . (19) 

With 7 = 5/3 for a monatomic gas, T oc (3 ~ 2 , and so the 
sound speed c a oc P 1 / 2 oc /3 _1 . The sound crossing time 
a/c s oc /3 2 . The fireball becomes almost immediately super¬ 
sonic. We can derive the behaviour of the sonic point from 
the definition that it is the position where the expansion 
rate matches the sound speed. Thus, 

1 _ v(r s ,t) _ Hr Si o _ Hr a _ Hr s 
C s (t) C a (t) /3c a (t) Cs,0 ' 

Assuming a uniform temperature profile and 7 = 5/3, this 
gives 


( 20 ) 


Cs,0 

r s = — = a 0 . 


( 21 ) 


Hence, we can see that the sonic point remains fixed in space 
at do. The initial density structure emerging through the 
sonic surface at radius do is subsequently ‘frozen in’ by the 
rapid expansion. 


2-4-2 Radiative Cooling Timescales 


The time required for the fireball to radiate away its in¬ 
ternal energy depends crucially on the temperature of the 
photosphere (assuming that the fireball is optically thick). 
For a uniform temperature distribution the timescale is very 
short: 


U_ 

U 


Mk 

m.HaR 2 T 3 


( M \ 

( R ) 

-2 r 

W 0 17 k g ; 

UO 8 mj 

U0 4 K ) 


( 22 ) 

(23) 


We would thus expect the region from which the radiation 
escapes to cool rapidly until the region becomes optically 
thin and the cooling efficiency drops. 

If the photospheric temperature T p is significantly 
cooler than the optically thick region, the radiative timescale 


U MkT 

U mn<jR 2 T p 



(24) 


(25) 


can become comparable to the evolutionary timescale of the 
fireball. 


2-4-3 Radiative Cooling Front 

Following the exact behaviour of radiative cooling and heat¬ 
ing is a complex process requiring us to follow the behaviour 
of the internal radiation field of the fireball through the 
optically thick to optically thin transition. Such detailed 
modelling is beyond the scope of the initial investigations 
presented here. We n ote that studies of cooling from op¬ 
tically thin plas mas (Lynden-Bcll & Tout 2001. Dalgarno 
fc McCray 1972| ) show that there is a significant reduction 
in efficiency due to hydrogen recombination and the con¬ 
sequent loss of free electrons below approximately 10 4 K. 
Since the photosphere is the region where the optical depth 
for escaping radiation is unity, this suggests that the photo¬ 
sphere may adopt a roughly constant effective temperature 
T p ~ 10 4 K. The relatively dense material in the fireball 
will be able to cool rapidly down to this temperature before 
the cooling rate stalls. We approximate a radiatively cooling 
fireball by maintaining a hot core region which cools adia- 
batically as it expands. Once gas elements cross the bound¬ 
ary to this region, we assume them to cool immediately to 
the effective temperature of the photosphere and thereafter 
adiabatically. We derive the radius of the photospheric sur¬ 
face below, using the blackbody luminosity and the thermal 
energy of the central region. 

To find the photospheric radius r p = ao/3p P in our ra¬ 
diative model we calculate the total thermal energy in a 
sphere of radius r p 


U = 


3 pkT 

2 pmn 


Anr 2 dr 


6-7T k ToaoPc.o 


prim /^O- 1 ) J 0 

Hence, differentiating, 

dU 6nkT 0 alp c 0 2 2 

V P e ,p 


rj~ dp. 


(26) 

(27) 


dp p /rrriH/? 3 ^ 7-1 ^ 

Radiative cooling at fixed effective temperature T p gives 


(28) 


dU . . . 2 rp, 2 rri4 

— = -4nHa 0 (3 y p aT p 
and so, using 


(29) 


dy P _ (dU 
d/3 ~ \ d/3 


dU 

dp P 


(30) 
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Photosphere Evolution 


Photosphere Evolution 



Figure 2. The evolution of the temperature boundary for the 
radiative model with M = 6.8 X 10 16 kg, ag = 9.6 X 10 7 m and 
To = 18 000 K. The solid line is for T p = 10 000 K and the 
dashed line for T p = 18 000 K. In each case we have initial values 
for rjp = 2, 3,4, 5. 


we have the differential equation for the core radius as a 
function of time 



We solve equation (|3^) numerically and plot the behaviour 
for typical parameters in Fig. [ 2 ]. We see how, in the La- 
grangian coordinate r/, the boundary migrates inward con¬ 
tinually. In Eulerian coordinates, plotted in Fig. [j, the 
boundary is initially advected outwards with the flow. When 
the inward migration exceeds the expansion rate, the pho¬ 
tosphere turns around and begins to collapse. In both plots 
the initial conditions for r/ p have relatively little impact on 
the behaviour at late times. Models with r) p <; 4 initially 
are indistinguishable as they rapidly converge on the same 
curve. Lower initial values for ri P also converge, albeit more 
slowly, on the same evolution. 

We plot the temperature profiles at different times for 
the radiative model in Fig. |dj. The profiles show how the ma¬ 
terial outside the core region cools very rapidly with distance 
away from the boundary. As a result, the shell of material 
with significant temperature outside the core is very thin. 
For computational convenience we set a minimum tempera¬ 
ture for any of the gas at 1 000 K. 

2.5 Ionization Structure 

In the LTE approximation, the density and temperature 
determine the ionization state of the gas at each point 
in space and time through the solution of a network of 
Saha equations. Atomic level populations are similarly deter¬ 
mined through Boltzmann factors and partition functions. 
For purely adiabatic cooling, the temperature retains its ini¬ 
tial uniform spatial profile, but decreases with time. The 
LTE ionization is therefore higher in the outer low-density 


Figure 3. The evolution of the temperature boundary in inertial 
coordinates for the radiative model with M = 6.8 X 10 16 kg, 
ao = 9.6 X 10 7 m and To = 18 000 K. The solid line is for 
T p = 10 000 K and the dashed line for T p = 18 000 K. In each 
case we have initial values for r/ p = 2, 3,4, 5. 


Temperature Profiles for Radiative Cooling Model 



V 


Figure 4. The temperature profile for a radiative fireball with 
M = 6.8 x 10 16 kg, a 0 = 9.6 x 10 7 m, T 0 = 18 000 K and 
T p = 18 000 K. Time runs from j3 = 0.1 for the hottest profile to 
/3 = 2.0 for the coolest in steps of’ 3 = 0.1 

regions which also move fastest. Thus in the LTE model 
high ionization emission lines are predicted to have broader 
velocity profiles. 

Once we have determined the evolution of T and p with 
time, we can follow the evolution of the ionization structure 
and determine the boundaries between different ionization 
states of various species. We do so for hydrogen and helium 
for an adiabatic Population II model in Fig. H and for carbon 
in Fig. j(| We plot boundaries for the same species for an 
isothermal Population II model in Figs. |t] and In both 
cases we use M = 10 18 kg, To = 18 000 K and ao = 10 s m. 

At a given time, the various ionic species form onion 
layers of increasing ionization on moving outwards from the 
fireball centre. The different temperature evolutions signifi¬ 
cantly alter the ionization evolution. The adiabatic temper¬ 
ature evolution causes the ionization states to alter rapidly 
throughout the structures when the fireball passes the ap- 
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Populotion II Adiabatic Ionization Boundaries 



Figure 5. Boundary between dominant species of hydrogen 
(solid) and helium (dashed) as a function of time for an adia¬ 
batic Population II model. 


Populotion II Adiabatic ionization Boundaries 



Figure 6. Boundary between dominant species of carbon as a 
function of time for an adiabatic Population II model. 


Population II Isothermal Ionization Boundaries 



Figure 7. Boundary between dominant species of helium as a 
function of time for an isothermal Population II model. Hydrogen 
is fully ionised throughout the calculation. 


Populotion II Isothermal Ionization Boundaries 



Figure 8. Boundary between dominant species of carbon as a 
function of time for an isothermal Population II model. 


propriate critical temperatures. The isothermal evolutions 
show the gentler dependance of ionization on density. The 
spatial structure remains roughly constant and only evolves 
slowly as the density drops. 


2.6 Ionization Freezing 

For the initial models presented here we have assumed that 
the ionization structure produced at any point is given by 
LTE; an assumption we justify in Appendix |X]. In reality, of 
course, once the LTE boundary is crossed, the ionization will 
drop as the net recombination rate gradually ekes away at 
the ions and we make a transition towards coronal equilib¬ 
rium. However, since we know the time-evolution of density 
and temperature for any region outside of the LTE limit¬ 
ing radius, we can derive the time and conditions when it 
crossed this LTE boundary. In principle then, we can follow 
the non-LTE evolution of the gas through this outer region 
and beyond the ion-electron equilibrium radius to the non¬ 
equilibrium ionization states present in the outer fringe. 

In a purely LTE fireball we expect high ionization states 
to occur in the outer regions as a result of the lower den¬ 
sity. With the more detailed prescription above, the higher 
ionization would result from the higher temperature of the 
fireball when the region in question crossed the LTE bound¬ 
ary. In practice, we expect that this correction may have 
little effect on the optical spectra because the emission is 
dominated by the higher density inner regions of the fire¬ 
ball. This effect may well become important, however, when 
we come to consider the ultraviolet spectra where high ion¬ 
ization lines are present; and which we hope to consider in 
a future paper. An alternative explanation for the ultravio¬ 
let behaviour may be that there is a hot, low density outer 
fringe to the fireballs. These conditions may occur in regions 
that have always been optically thin and inefficient at cool¬ 
ing radiatively. These lines may then arise from a transition 
region similar to that between the solar chromosphere and 
corona. 
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3 RADIATIVE TRANSFER 

The radiative transfer equation has a formal solution 

1 = J Se~ T dr, (33) 

where I is the intensity of the emerging radiation, S is the 
source function, and r is the optical depth measured along 
the line of sight from the observer. This integral sums contri¬ 
butions Sdr to the radiation intensity, attenuating each by 
the factor e~ r because it has to pass through optical depth 
r to reach the observer. 

Since we have assumed LTE, the source function is the 
Planck function, and opacities both for lines and continuum 
are also known once the velocity, temperature, and density 
profiles and element abundances are specified. The integral 
can therefore be evaluated numerically, either in this form 
or more quickly by using Sobolev resonant surface approxi¬ 
mations. 

The above line integral gives the intensity I(y) for lines 
of sight with different impact parameters y. We let x mea¬ 
sure distance from the fireball centre toward the observer, 
and y the distance perpendicular to the line of sight. The 
fireball flux, obtained by summing intensities weighted by 
the solid angles of annuli on the sky, is then 

fW = j o I{\y)^dy (34) 

where d is the source distance. 


3.1 Continuum Lightcurves and Spectra 


In a forthcoming paper (Pearson, Horne & Skidmore, in 
prep.) we will consider generic analytic models for fireball 
behaviour applicable to a variety of s ystems. We use the 
general opacity given in equation (A4). We quote here an 
approximation to the lightcurve behaviour being the sum of 
optically thick and optically thin contributions: 


fu = 

P,oPo t oS (TO) 


(35) 

where 
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Here /3g is the time that the fireball first becomes optically 
thin through its centre and to is the current optical depth 
through the fireball centre. 

The peak flux occurs when 


To 


(42) 


Predicted Adiabatic Lightcurves 



0 0.2 0.4 0.6 0.8 1 


P 


Figure 9. The predicted fluxes at 3 000 (lowest), 4 000, 5 000, 
6 000, 7 000 and 8 000 A (highest) for an adiabatic fireball 
with M = 10 17 kg and To = 10 000 K. The upper family have 
ao = 10 8 m and the lower ao = 3 X 10 7 m. The dashed lines show 
the rapidly declining optically thick contribution at 8 000 A. 


which occurs at a time 

/3pk = /3 0 e-A (43) 

We note that equations (37) and (^) predict that the peak 
flux will occur at different times for different wavelengths. 
In principle we could use this as a test for the effective I\ 
However, in practice, both the assumptions in our treatment 
and the necessary accuracy of observations, might make it 

3 

difficult to detect the difference between /3 P k ~ As for an 
isothermal fireball and /3 p k ~ A* for adiabatic cooling. In 
contrast, the late time behaviour is predicted to have a more 
significant dependance on T 

( 44 ) 

Coverage of a full flare may well make the /? ~ f3~ 2 or /3~ 3 
difference discernable. 

The total flux at several wavelengths and the optically 
thick contribution at 8 000 A are plotted in Figs. H and jlOj 
for 2 sets of physical conditions, using Population I compo¬ 
sition and a distance of 100 pc. The corresponding spectral 
distributions for one of the parameter sets are plotted in 
Figs. 0 and 0 We would expect the early and late time 
behaviours to be accurately represented by this analysis, but 
the detailed behaviour when both free-free and bound-free 
opacity sources are important is less exact. 

The lightcurves for the isothermal models show a rise to 
the peak value before turning over and undergoing a slower 
decline. The predicted fluxes are very similar to observed 
fireballs and morphologically the similarity to observations is 
remarkable. It is also clear that the peak flux is independent 
of our choice of fiducial size. We can see how the combination 
of the overall spectral slope and the large Balmer jump at 
early times leads to the 3 000 A curve exceeding that for 
4 000 A. As ionization occurs due to the decreasing density, 
bound-free opacity becomes less important and the Balmer 
jump declines. The 3 000 A curve then drops below the 
4 000 A and behaves similarly to the other curves in the 
family. The size of our Balmer jump, however, is still smaller 
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Predicted Isothermal Lightcurves 



0 0.5 1 1.5 2 


Figure 10. The predicted fluxes at 3 000 (lowest), 4 000, 5 000, 
6 000, 7 000 and 8 000 A (highest) for an isothermal fireball with 
M = 10 17 kg and To = 10 000 K. The family on the left have 
ao = 10 8 m and the right ao = 3 X 10' m. The dashed lines show 
the optically thick contributions at 8 000 A 


Predicted Adiabatic Spectral Distributions 



Wavelength (A) 


Figure 11. The predicted spectral distribution for an adiabatic 
fireball with M = 10 17 kg, To = 10 000 K and ao = 3 X 10' m. 


Predicted Isothermal Spectral Distributions 



Wavelength (A) 


Figure 12. The predicted spectral distribution for an isothermal 
fireball with M = 10 17 kg, To = 10 000 K and ao = 3 X 10 7 m. 


than that observed suggesting that we may need to improve 
the treatment of opacity across this region. 

The adiabatic lightcurves, however, show only a decline 
from a large initial flux. We would expect the evolution of 
the optically thick region here to be similar to that of the 
radiative core photosphere ie. it would be initially advected 
outward with the flow before eventually turning over and 
collapsing back to zero size. The emitted flux, however, does 
not follow this scheme as the fireball is extremely hot at early 
times and the evolution of temperature and the blackbody 
intensity dominates the lightcurve evolution. The spectra in 
this case only begin to show a Balmer jump late on in the 
evolution once the fireball has cooled enough for significant 
recombination to have taken place. 

3.2 Fireball Tomography 

For our spherical expansion the component of velocity along 
the line of sight is proportional to the distance from the 
centre towards the observer, ie. Vx = HX. As a result, 
each velocity in the emission line profile originates from a 
unique surface perpendicular to the line of sight. This is 
illustrated in Fig. 0 which shows a side view of the fireball 
with the contribution to the flux at 3 different velocities in 
the Ha profile. On the blue side of line centre, the resonant 
surface, approximately perpendicular to the line of sight, 
is on the near side of the fireball. The fireball is optically 
thick through its centre but sufficiently transparent at larger 
impact parameters so that small continuum contributions to 
the flux are made from a peanut shaped region along the y- 
axis about the origin. 

The resonant surface at the line centre passes through 
the centre of the fireball. There is a continuum contribution 
from the rising density in front of the resonant plane around 
the origin and then a line contribution from the resonant 
position. The larger optical depth through the centre com¬ 
pared to larger impact parameters gives rise to the curved 
shape. The contribution from the back of the resonant region 
through the line centre is unable to escape to the observer 
through the overlying material giving rise to the distorted 
shape. 

The flux on the red side of the line shows the two con¬ 
tributing regions more clearly. A continuum contribution 
with the density peak and then another peak at the res¬ 
onant position. 

Such fireball tomography may allow us, in future, to 
probe the structure of a fireball using high-resolution spec¬ 
tra. For example, in Fig. [b|, the Gaussian velocity profile 
arises from the assumed Gaussian density profile in equa¬ 
tion ©■ A different density profile would change the shape 
of the line profile. 

3.3 Numerical Considerations 

Our numerical integrations of equation ( ^ ) were carried out 
at appropriate wavelength intervals to resolve the emission 
line profiles and with spatial steps to reflect the variations 
in the density and ionization structures. For conservative 
parameters T = 10 4 K and Ao = 3 500 A, the thermal line 

width is AA ~ ^ {-ff) 2 « 0.1 A. Opacity information was 
therefore generated and stored at wavelength and spatial 
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Figure 13. Contribution to the flux at the wavelengths indicated from different positions in the fireball. The observer is at +oo 
on the X axis and Y is the impact parameter. Both axes are labelled in multiples of the current lengthscale a = 4.2 X 10 7 m. 


intervals of 0.03 A and 0.1 a respectively. When generating 
the full spectra, intensities were calculated every 1 A, giving 
a velocity resolution of 86 km s _1 at 3 500 A. The integration 
was carried out with spatial steps that were calculated to be 
either At = 0.005 or Air = 0.1 a whichever was shorter. 

In general the opacity at any position in the fireball 
consisted of both continuum and line contributions. Contin- 
uu m opacities were generally calculated using the me thods 
of Gray (1976 ) with the exception of bound -free (Celh 

man 19 52 ) and free-free ([g tilley fe Callawa y 197C ), He~ ( [Mc¬ 
Dowell. Williamson fe Myerscough 1966) and Hel bound-free 


(Huang 1948). Line opacities used oscillator strengths and 
energy level information downloaded from the NIST atomic 
database!!. 

We aim, in the longer term, to speed-up our code by 
using a Sobolev method, ie. use the fact that velocity is a 
monotonically increasing function of position, to find the 
resonant point in the fireball at each wavelength for a given 
line. We might also achieve time savings by identifying re¬ 
gions where no lines make significant contributions to the 
flux. These methods would enable us to use coarser grids 
or precalculated tables when appropriate. This would sim¬ 
plify the integration scheme and possibly allow us to fit both 
the emission lines and continuum simultaneously with an 
amoeba code. 


4 COMPARISON WITH OBSERVATIONS 


4.1 Observed Flare in AE Aqr 

To test our fireball m odel, we aimed to repr oduce the ob¬ 
servations reported in Skidmore et al. (2002) of high time- 
resolution optical spectra of an AE Aqr flare taken with the 
Keck telescope. The lightcurve formed from these data is 
plotted in Fig. |lj. Although the dataset only lasts around 
13 minutes it shows the end of one flare and beginning of 
another. The major features of the flaring observed in the 
system are apparent in the lightcurve. 

We extracted the spectrum of the second flare by sub¬ 
tracting the quiescent spectrum (Q) from the spectrum at 
the top of the new flare (P). The resultant optical spectrum 
is plotted in the bottom panel of Fig. 0 Gaussian profiles 
fitted to several lines gave the parameters in Table As¬ 
suming the instrumental profile is also Gaussian, we remove 
the formal instrumental blurring (9.8 A FWHM) to estimate 
the expansion velocity. 

These widths suggest a velocity dispersion of around 
700 km s _1 . Ascribing this to the velocity at la of the den- 

at ao. Such 

-l 


sity distribution gives a velocity of ~ 950 km s 
a high velocity is larger tha n the expected V ~ 300 km s 
from Doppler tomography (Welsh, Horne & Corner 199£) 


§ http://physics.nist.gov/cgi-bin/AtData/main_asc 


and may indicate that an additional, unknown source of 
blurring is present. Examination of the arc line spectra sup¬ 
port this suggestion. Since these were taken by flooding a 
large aperture with light it is not possible to accurately mea¬ 
sure the instrumental effect from these lines, but the line 
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Cornpcmson of Voce one Dole Spec:ro 



3 <2 


Observed Continuum Lightcurve 



12.9 12.95 13 13.05 13.1 

Time (UT) 


Figur e 14. Observed lightcurve reported in Skidmore et al. 


(2002) ) 


showing the decline of an old flare and start of a small 


new flare. 


edges do suggest that the instrumental blurring is signifi¬ 
cantly larger than the formally derived value. 


4.2 Spectral Fit at Peak of Flare 

To estimate the fireball parameters M, T and a from the 
observed optical spectra we began with consideration of the 
continuum flux. We considered both Population I and Pop¬ 
ulation II abundances. Population I abun dances were se t 
using the solar abundances presented by Dappen (2000 ). 
For population II composition, we decreased from 0.085 to 


Line 

Fitted FWHM 
(km s -1 ) 

Deconvolved FWHM 
(km s -1 ) 

Ha 

1909 ± 7 

1856 ± 7 

Hel 5876 

1510 ±6 

1424 ± 6 

Hel 4922 

1353 ± 13 

1214 ±13 

H/3 

1681 ± 3 

1569 ± 3 

Hy 

1743 ± 3 

1589 ± 3 

H5 

2171 ±4 

2050 ± 4 

He 

1645 ± 5 

1469 ± 5 

CaH 

2423 ± 21 

2306 ± 21 


Table 2. Fitted and deconvolved widths of several emission lines. 


0.028 the ratio of the helium to hydrogen number densities 
and r educed the ratio for each metal abundance to 5% of 
solar ( Bowers fe Deeming 1984 ). In either case we used an 
amoeba algorithm (Press et al. 1986) to optimize the fit by 
reducing \ 2 to a minimum for the 6 measured continuum 
fluxes in Table The best fit parameters M, T and a and 
related secondary quantities are listed in Table 

Using Population I composition, we derived parame¬ 
ters for the fireball at its peak of M — 3.7 x 10 16 kg, 
do = 5.1 x 10 7 m and To = 17 000 K. Such a mass rep¬ 
resents approximately 100 s of the Li mass transfer rate 
derived from equation (:2j). Fig. [T^ and Table [l] show that 
these parameters provide a reasonably good fit to the ob¬ 
served Paschen continuum and Balmer Jump. 

Population II composition produced best fit parameters 
of M = 6.8 x 10 16 kg, ao = 9.6 x 10 7 m, T 0 = 18 000 K. 

Both compositions arrive at a model fit with remark¬ 
ably similar calculated fluxes which have a slightly steeper 
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Quantity 

Pop I 

Pop II 

M (10 16 kg) 

3.7 

6.8 

ao (10 7 m) 

5.1 

9.6 

To (K) 

17 000 

18 000 

u(ao) (km s 1 ) 

170 

170 

l/H (s) 

300 

560 

Po (10 8 kg m 3 ) 

5.1 

1.4 


Table 4. Parameters used for the model spectra and useful sec¬ 
ondary quantities. 


slope to the Paschen continuum than the observations. The 
fits also produce a large Balmer jump but appears to have 
difficulty in creating one quite as strong as that observed. A 
change in the 3rd significant figure of one of the fit param¬ 
eters generally leads to a similar change in the calculated 
fluxes. Both sets of parameters have similar temperatures 
but the Population II models require roughly a factor 2 in¬ 
crease in mass and lengthscale equivalent to a factor 4 de¬ 
crease in density. The central densities implied by both sets 
of parameters are consistent with the typical values for the 
gas stream derived in section 2.1 and with the temperature 
range and emitting area of Beskrovnaya et al.( 1996). 

Fig. § shows the area around the H/3 line produced 
with the Population I best fit parameters and with zero and 
modest expansion velocities, showing the intrinsic thermal 
resolution of the lines and the Doppler blurring due to the 
expansion before instrumental effects are added. Even with 
a relatively modest fiducial expansion velocity of 100 km s -1 
we can see that many of the weaker lines are blurred into 
obscurity simply by the Doppler effect. The H/3 line at 
4 861 A appears to be optically thick since its peak flux 
decreases only slightly as the line width increases. In con¬ 
trast, the Aril line at 4 807 A decreases greatly with Doppler 
broadening and forms a blend with the weaker blueward Nil 
line indicating that it is optically thin. 


tmg 1 

In Fig. llT] we plot the observed data around H/3 along¬ 
side Population I models generated with various expansion 
velocities. We consider the H/3 region because the observed 
Ha line shows a complex structure. We blurred the model 
spectra with a Ga ussian of 9.8 A FWHM and assumed a 
distance of 102 pc (Friedjung 1997). To match the peak flux 
of H/3 requires an expansion velocity of 75 km s -1 at do. 
The Hel line at 4 923 A suggests an expansion velocity 
of 40 km s _1 . The model lines are clearly narrower than 
the data. Considering this and the possibility of additional, 
unidentified blurring suggests that a more reliable estimate 
might be the integrated line flux. In this case, the H/3 and 
Hel line suggest expansion velocities of 240 km s^ 1 and 
100 kin s -1 respectively at a o. Consequently, we adopt an 
expansion velocity at ao of 170 km s _1 for the simulations. 

The optical spectra produced with both compositions 
are compared to the observed mean spectrum at the flare 
peak in Fig. [1^. We see that the parameters, derived from 
the peak continuum fluxes alone, allow us to derive opti¬ 
cal spectra with integrated line fluxes and ratios compa¬ 
rable to the observations. Both models have line widths 
narrower than the observations, and hence peak line fluxes 
greater than those observed, which may result from ad¬ 
ditional unidentified instrumental blurring. The Popula- 


Comparison of Data to Various Model Expansion Velocities 



Figure IT. Observed data around H/3 plotted alongside model 
spectra with Population I composition and expansion velocities 
at ao running from 0 to 110 km s —1 at 10 km s —1 intervals. 


tion I models show the fiat Balmer decrement of satu¬ 
rated Balmer lines apparent in the data. Comparing the 
lines that are present, however, particularly in the 4000 — 
5000 A range, the Population II model appears to pro¬ 
duce a much better fit: suppressing the metal lines which 
do not appear in the observed spectra. Specifically, we note 
that SIIIIA3080, CIIA4268, the SiIIIAA4554, 4555, 4569, 4576 
blend, OIIA4650, the SiIA5667, SiIIIAA5697, 5698, 5741 and 
ArIIA5726 blend and CIIAA7233, 7238 are present in the 
Population I models but are not present or are far weaker 
in the data. Only the weak CIIA3922, CIIA5147, the 
SiIIAA6349, 6373 blend and the OIA7774, 7775, 7778 blend 
are not reproduced by the Population II results. 


4.3 Lightcurves and Spectral Evolution 

The timescale for the evolution of the fireball implied by the 
observed lightcurve can be reproduced by the parameters 
derived solely from fitting to the peak spectrum. The evo¬ 
lutionary timescale ao/v(ao) = 1 /H is approximately 300 s 
and 560 s for Population I and II parameters respectively. 
The encouraging result that these values are similar to the 
observed timescale suggests that our fireball models are on 
the right track. 

We thus consider the time evolution of the fireball mod¬ 
els in more detail. Figs, [is] and ^9[ for Population I and Pop¬ 
ulation II abundances respectively, present the lightcurves 
and spectral evolution derived under the three assumptions 
we consider for the thermal history of the fireball: adiabatic, 
isothermal, and radiative. In each case, the model fireball is 
constrained to evolve through the state we found from the 
fit to the observations at the peak of the flare. This moment 
is assigned to t = 0 and hence our assignment of 0 subscripts 
for the fitted parameters above. 

The isothermal models produces a lightcurve that rises 
and falls, as observed. Moreover, the peak in this lightcurve 
is in fact reasonably close to the time when the spectrum 
matches the spectrum observed at the peak of the flare. 
The isothermal fireball spectra, both on the rise and on the 
fall of the flare, exhibit strong Balmer emission lines and a 
Balmer jump in emission. Thus the isothermal fireball model 
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Range 

(A) 

Central Wavelength 

(A) 

Mean Flux 
(mjy) 

Pop I Fluxes 
(mJy) 

Pop II Fluxes 

(mJy) 

3602-3630 

3616 

30.68 ± 0.20 

27.3 

28.3 

4150-4250 

4200 

6.812 ±0.019 

6.65 

6.64 

4520-4640 

4580 

7.442 ±0.017 

7.52 

7.49 

4735-4780 

4758 

7.394 ± 0.028 

7.89 

7.85 

5400-5600 

5500 

9.114 ±0.016 

9.13 

9.13 

5600-5775 

5688 

9.346 ±0.018 

9.37 

9.39 


Table 3. Mean flux measured in the 6 continuum regions with the calculated continuum fluxes from the models. 


Examination of Expansion Blurring 


o 



Wavelength (nm) 


Figure 16. Model spectra around H/3 for Population I composition with zero expansion velocity (solid) and an expansion velocity 
of 100 km s -1 at a o (dotted) before instrumental and distance effects are added. 


provides a good fit not only to the peak spectrum but also 
to the time evolution of the AE Aqr flare. 

The adiabatic and radiative models both fail to repro¬ 
duce the observed lightcurve or spectra. The adiabatic model 
matches the observed spectrum at one time, but is too hot 
at early times and too cool at late times. As a result the 
lightcurve declines monotonically from early times rather 
than rising to a peak and then falling, and the spectra have 
lines that do not match the observed spectra. This disap¬ 
pointing performance is perhaps not too surprising, since 
the temperature in this model falls through a wide range 
while the observations suggest that the temperature is al¬ 
ways around 10 4 K. 

For the adiabatic fireball, equation (^) implies that 


the expansion rate is a measure of the initial temperature. 
Hence, for an adiabatic fireball we can derive its age from the 
time taken to cool to the fitted temperature. For the Popu¬ 
lation I parameters, 170 km s _1 gives an initial temperature 


of 1.25 x 10° K, an age of 260 s and hence from (|l 
tial lengthscale of 5.9 x 10 6 m. The Population II parameters 
give an age of 494 s and initial lengthscale of 1.1 x 10' m. We 
note, however, that the lightcurve for the adiabatic fireballs 
do not show the rise to a pe ak for the same reason as the 
theoretical curves in section 3.1 To explain the rise phase 
of the observed lightcurve we must interpret it as resulting 
from the collision of the initial gas blobs. Under this assump¬ 
tion, we can estimate the age and initial size independently , 
for the adiabitic evolution, directly from the data; knowing 
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Comparison of Popu'cfon Mode Spectra one Continuum ,'gntcjrves 



3000 4000 5000 6000 7000 8000 -200 -100 0 100 200 300 

Wavelength (A) Time (s) 


Figure 18. Population I spectra at t = —239, —179, —119, —60,0 s, (/3 = 0.2, 0.4,0.6, 0.8,1.0) for the three cooling laws alongside 
their respective continuum lightcurves at 5 350 A. 


Compa^son o' Pooo : ci"on l Voa^ Specuo and Continuum jghtc,. r ves 



3000 4000 5000 6000 7000 8000 -400 -200 0 200 400 600 

Wavelength (X) Time (s) 


Figure 19. Population II spectra at t = —450, —337, —225, —113, 0 s, (/3 = 0.2, 0.4, 0.6,0.8,1.0) for the three cooling laws alongside 
their respective continuum lightcurves at 5 350 A. 


© 0000 RAS, MNRAS 000, 000-000 









Fireball Models for Flares in AE Aquarii 15 


the rise time of the flare and the typical closing velocities: 

tage trise — 63 S (45) 

l ~ Vt liae = 3 x 10 5 x 60 = 1.8 x 10 7 m. (46) 

The radiative model looks roughly the same as the adia¬ 
batic model. This occurs because the cooling front (as shown 
in Fig. [i]) occurs well outside the photosphere. We therefore 
“see” into the adiabatically-cooling core. This is not con¬ 
sistent with our assumption that the rapid cooling occurs 
near the photosphere as a result of the opacity drop there 
and suggests that a more sophisticated approach involving 
heating of the gas by emergent photons may be required for 
an accurate treatment of such a model. 


5 DISCUSSION 


The spectra produced with Population II composition re¬ 
produce the optical observations well. The fitted parameters 
are consistent with both the expected conditions in the mass 
transfer stream and with the lower limits on the density im¬ 
plied by the uv observations. The presence of high ionization 
and semi-forbidden lines in the uv data can be understood 
qualitatively in terms of their formation in the low density 
outer fringes of the fireball and the wide variety of ioniza¬ 
tion states by the wide range of densities in the expanding 
fireball structure. For f3 = 1, /x = 0.53 and Population II 
parameters, the density of 10 17 m -3 , important for the uv 
semi-forbidden lines, occurs at r/ ~ 2.3. From Fig. A1 we see 


that this lies outside the limit for LTE behaviour. Quantita¬ 
tive fits for the fringe and uv behaviour, therefore, remain to 
be addressed in the non-LTE regime extension to this work. 

Parameters implying a lower overall density would re¬ 
sult in an optically thinner fireball and spectrum. The 
Balmer decrement would increase as would the equivalent 
widths of the optically thick Balmer lines as the continuum 
dropped away from the blackbody envelope. Continued de¬ 
crease would result in the lines eventually becoming optically 
thin and losing the characteristically consistent strengths 
present when they are all saturated. 

The behaviour of the lightcurve for an isothermal 
model, like its theoretical counterpart is much more simi¬ 
lar to the observed lightcurve than the other models giving 
a peak close to the observed time without the need to invoke 
the collision process. Clearly though, an expanding gas ball 
would be expected to cool both from radiative and adiabatic 
expansion effects. 

The radiative model shows very similar behaviour to 
the adiabatic models. The thin shell of cooling material at 
around 10 4 K does not provide sufficient opacity to mimic 
the isothermal behaviour and the adiabatic core is still visi¬ 
ble to the outside observer. This effect appears to be insensi¬ 
tive to the exact choice of temperature we assign to material 
emerging from the central region. 

We noted in section 2.6 that an accurate picture of the 
ionization would need to follow the LTE to non-LTE ion¬ 
ization transition. However, even considering only the early 
evolution of the fireballs, the isothermal lightcurve and spec¬ 
tral behaviour is stiller closer to the observations than the 
adiabatic. 

In short, we find that isothermal fireballs reproduce the 


observations rather well, whereas adiabatically cooling fire¬ 
balls fail miserably! How can the expanding gas ball present 
a nearly constant temperature at its surface? 


5.1 Thermostat 

We offer two possible mechanisms which may operate to 
maintain the apparent fireball temperature in the 1 -2 x 
10 4 K region. The first relies on the fact that both free-free 
and bound-free opacity decrease with temperature and thus 
opacity peaks just above the temperature at which hydro¬ 
gen begins to recombine. In the model where a core region 
is cooling through radiation from its surface, we can envi¬ 
sion a situation where the material just outside the core has 
reached this high opacity regime and is absorbing significant 
amounts of energy from the core’s radiation field. This would 
help to counteract the effect of adiabatic cooling and may 
maintain an effectively isothermal blanket around the core. 
If the temperature in the blanket were to rise, the opacity 
would decrease, more energy would escape and the blan¬ 
ket would cool again. Similarly, until recombination, cooling 
would result in higher opacity and therefore greater heating 
of the blanket region. For such a thermostatic mechanism 
to work, the optical depth through the blanket region would 
need to be high which would be consistent with it being the 
photosphere as seen by an outside observer. Realistic simu¬ 
lations of such a model would require a detailed treatment 
of the radiation field and its heating effect using a more 
sophisticated method than the one we have employed here. 


5.2 External Photoionization 

Alternatively the ionization of the fireball might be held at 
a temperature typical of the 1-2 x 10 4 K range through 
photoionization by the white dwarf. In Fig. f2(] we plot 
the hydrogen ionization and recombination rates as a func¬ 
tion of temperature for a typical electron number density 
tie = 3 x 10 18 m -3 using additional routines for photoion¬ 
ization from Verner et al. (1996) and collisional ionization 
by Verner using the results of Voronov (1997). We calculate 
a conservative overestimate for the self-photoionization of 
the fireball assuming the sky is half filled by blackbody ra¬ 
diation at the given fireball temperature. Also shown is the 
photoionization due to a blackbody with the same radius 
and distance as the white dwarf. The general white dwarf 
temperature is uncertain but, as mentioned earlier, there is 
evidence to suggest that the hot spots have a temperature 
of around 3 x 10 4 K. We can see that the flux from the white 
dwarf overtakes that from the fireball itself as the dominant 
ionization mechanism at ~ 10 4 K. We have assumed a typ¬ 
ical distance from the white dwarf for the fireball equal to 
that of the LI point from the white dwarf (although in a 
different direction). The recombination rates for this plot 
have been calculated assuming n e is a constant and so for 
temperatures below about 10 4 K are conservative overesti¬ 
mates of the LTE recombination rate that will occur in our 
fireball. In spite of this, the white dwarf photoionization at 
3 x 10 4 K comfortably exceeds the recombination rates down 
to at least 10 3 K and, hence, we would expect the fireball 
to remain almost completely ionized in this case. This sug¬ 
gests that the white dwarf radiation field may be the cause 
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Figure 20. Ionization and recombination rates for a hydrogen 
atom or ion respectively with n e = 3 X 10 18 rn "" :i . 

of the fireball appearing isothermal when we consider the 
lightcurve behaviour and consistent line strengths and ra¬ 
tios. 

Ionization by an external black body at a different tem¬ 
perature will clearly cause our fireball to no longer be in 
LTE and require a non-LTE model to follow the fireball 
behaviour. However, we can anticipate the results of a full 
non-LTE treatment as follows. Writing the total ionization 
rate per atom per second from all processes as b(nu,n e ,T) 
and similarly the total recombination rate as a(n e ,T), we 
have the simple equilibrium condition 

bn n = an i. (47) 

Hence the ionization fraction is given by 


- — 1 

n n + rii b + a 

We calculate the evolutions of the boundary between 
HI and HII dominated regions for a pure hydrogen fireball 
using M = 10 16 kg, ao = 10 7 m and To = 18 000 K. Since in 
this simplified case n e = nr, we can easily iterate to produce 
self-consistent ionization profiles in cases where b contains 
contributions from ordinary atomic processes and also an 
additional contribution from white dwarf photoionization. 
We plot the evolutions in Fig. ^1] for no white dwarf pho¬ 
toionization and with rates appropriate to photoionization 
from a white dwarf at 2 x 10 4 K and 3 x 10 4 K. We see that 
the white dwarf can have a significant effect on the ioniza¬ 
tion structure of the fireball and keep large parts of it in an 
ionized state. 


6 SUMMARY 

We have shown for AE Aqr how the observed flare spectrum 
and evolution is reproducible with an isothermal fireball at 
T = 18 000 K with Population II abundances but not when 
adiabatic cooling is incorporated. We suspect that the cause 
of the apparently isothermal nature is a combination of two 
mechanisms. First, a nearly isothermal photosphere which 
is self-regulated by the temperature dependance of the con¬ 
tinuum opacity and the hydrogen recombination front and 


Irradiated Adiabatic HI/HII Boundary Radius v Time 



P 

Figure 21. Comparison of the evolution of the hydrogen ioniza¬ 
tion boundary for purely LTE and LTE plus white dwarf pho- 
toionization cases. 


second, particularly at late times, by the ionizing effect of 
the white dwarf radiation field. 

The observed photosphere in the isothermal models is 
initially advected outward with the flow before the decreas¬ 
ing density causes the opacity to drop and the photosphere 
to shrink to zero size. This gives rise to a lightcurve that 
rises and then falls. Emission lines and edges arise because 
the photospheric radius is larger at wavelengths with higher 
opacity. Surfaces of constant Doppler shift are perpendicular 
to the line of sight and thus Gaussian density profiles give 
rise to Gaussian velocity profiles. 

Direct application of our LTE fitting method to the uv 
region would not reproduce the observed spectra with their 
variety of ionization states and semi-forbidden and permit¬ 
ted lines. We have outlined techniques and possible improve¬ 
ments to the model that will increase our understanding 
further and enable us to take the work into the uv. More 
detailed modelling of the ionization and recombination pro¬ 
cesses is suggested that also incorporates the heating effects 
of both the white dwarf’s and the fireballs own photon fields. 
The self-consistent solution to these rate equations will ex¬ 
tend the treatment to the non-LTE regimes discussed in Sec¬ 
tion and Appendix |a| These modifications may in turn 
have consequences for the predicted optical spectra. 

Improved modelling of these fireballs offers us the 
chance to probe the chemical composition of the secondary 
star in AE Aqr. Normally this is difficult because the spec¬ 
trum of the relatively dim secondary star is contaminated 
by light from other components in the system. We expect 
that the fireball models will be applicable more generally to 
the flickering observed in most CV systems. 


Acknowledgements 

We would like to thank Gary Ferland and Kirk Korista for 
informative discussions regarding atomic ionization and re¬ 
combination processes and the routines for calculating them. 
We are grateful to the anonymous referee and editor for help¬ 
ful comments that improved the presentation of this paper. 


© 0000 RAS, MNRAS 000, 000-000 











Fireball Models for Flares in AE Aquarii 17 


REFERENCES 

Abada-Simon, M., Bastian, T. S., Horne, K. D., Robinson, E. L., 
Bookbinder, J. A., 1995, in Buckley, D. A. H., Warner, B., 
eds., Proc. Cape Workshop on Magnetic Cataclysmic Vari¬ 
ables, ASP Conf. Series, 355 
Arnaud, M., Raymond, J., 1992, ApJ, 398, 394 
Bak, P., 1996, How nature works: the science of self-organized 
criticality, Copernicus, New York 
Bastian, T. S., Dulk, G. A., Chanmugam, G. 1988, ApJ, 324, 431 
Bastian, T. S., Dulk, G. A., Chanmugam, G. 1988, ApJ, 330, 518 
Beskrovnaya N., Ikhsanov N., Bruch A., Shakhovskoy N., 1996, 
A&A, 307, 840 

Bowers, R. L., Deeming, T., 1984, Astrophysics I Stars, Jones and 
Bartlett, Boston 
Bruch A., 1991, A&A, 251, 59 
Bruch, A., 1992, A&A, 266, 237 
Bruch A., Gruuter M., 1997, Act Ast, 47, 307 
Cota, S. A., 1987, Ph.D. Thesis, Ohio State Univ. 

Dalgarno, A., McCray, R. A., 1972, ARA&A, 10, 375 
Dappen, W., 2000, in Cox, A. N., Allen’s Astrophysical Quanti¬ 
ties, Chapter 3, Springer-Verlag, New York 
de Jager, O. C., Meintjes, P. J., O’Donoghue, D., Robinson, E. 

L. , 1994, MNRAS, 267, 577 

Elsworth, Y. P., James, J. F., 1982, MNRAS, 198, 889 
Eracleous, M., Horne, K. D., Robinson, E. L., Zhang, E.-H., 
Marsh, T. R., Wood, J. H., 1994, ApJ, 433, 313 
Eracleous, M., Horne, K. D., 1996, ApJ, 471, 427 
Frank, J., King, A. R., Raine, D. J., 1992, Accretion Power in 
Astrophysics, Cambridge Univ. Press 
Friedjung, M., 1997, NewA, 2, 319 
Geltman, S., 1962, ApJ, 136, 935 
Glasco, H. P., Zirin, H., 1964, ApJS, 9, 193 
Goldston, R. J., Rutherford, P. H., 1995, Introduction to Plasma 
Physics, Bristol, IOP Publishing 
Gray, D. F., 1976, Observations and Analysis of Stellar Photo¬ 
spheres, New York, Wiley 

Horne, K. D., 1999, in Hellier, K., Mukai, K., eds, Annapolis 
Workshop on Magnetic Cataclymic Variables, ASP Conf. Se¬ 
ries, 157, 357 

Huang, S., 1948, ApJ, 108, 354 

Hameury, J. M., King, A. R., Lasota, J.-P., 1986, A&A, 162, 71 
Jameson, R. F., King, A. R., Sherrington, M. R., 1980, MNRAS, 
191, 559 

Kingdon, J. B., Ferland, G. H., 1996, ApJS, 106, 205 
Kovetz, A., Prialnik, D., Shara, M. M., 1988, ApJ, 325, 828 
Kuijpers, J., Fletcher, L., Abada-Simon, M., Horne, K. D., Raadu, 

M. A., Ramsay, G., Steeghs, D., 1997, A&A, 322, 242 
Landini, M., Monsignori Fossi, B. C., 1990, A&AS, 82, 229 
Landini, M., Monsignori Fossi, B. C., 1991, A&AS, 91, 183 
Lynden-Bell, D., Tout, C., 2001, ApJ, 558, 1 

McDowell, M. R. C., Williamson, J. H., Myerscough, V. P., 1966, 
ApJ, 144, 827 

Mazzotta, P., Mazzitelli, G., Colafrancesco, S., Vittorio, N., 1998, 
A&AS, 133, 403 

Meyer, F., Meyer-Hofmeister, E., 1983, A&A, 121, 29 
Papaloizu, J. C. B, Bath, G. T., 1975, MNRAS, 172, 339 
Patterson J., 1979, ApJ, 234, 978 

Patterson, J., Branch, D., Chincarini, G., Robinson, E. L., 1980, 
ApJ, 240, L133 

Pequignot, D., Petitjean, P., Boisson, C., 1991, A&A, 251, 680 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., Flannery, B. P., 
1986, Numerical Recipes in Fortran, Cambridge Univ. Press 
Ritter, H., 1988, A&A, 202, 93 
Sarna, M. J., 1990, A&A, 239, 163 
Shull, J. M., Van Steenberg, M., 1982, ApJS, 48,95 
Skidmore, W., O’Brien, K., Horne, K. D., Gomer, R., Oke, J. B., 
Pearson, K. J., 2002, MNRAS, submitted 


Stilley, J. L., Callaway, J., 1970, ApJ, 160, 245 

Thorstensen, J. R., Ringwald, F. A., Wada, R. A., Schmidt, G. 

A., Norsworthy, J. E., 1991, AJ, 102, 272 
van Paradijs J., Kraakman H., van Amerongen S., 1989, A&AS, 
79, 205 

Verner, D. A., Ferland, G. J., 1996, ApJS, 103, 467 
Verner, D. A., Ferland, G. J., Korista, K. T., Yakovlev, D. G., 
1996, ApJ, 465, 487 

Voronov, G. S., 1997, Atomic Data and Nuclear Data Tables, 65, 
1 

Warner, B., 1995, Cataclysmic Variable Stars, Cambridge Univ. 
Press 

Welsh W.F., Horne K., Oke B., 1993, ApJ, 406, 229 
Welsh, W. F., Horne, K. D., Gomer, R., 1995, MNRAS, 275, 649 
Welsh, W. F., Horne, K. D., Gomer, R., 1998, MNRAS, 298, 285 
Wynn, G. A., King, A. R., Horne, K. D., 1995, in Buckley, D. 
A. H., Warner, B., eds., Cape Workshop on Magnetic Cat¬ 
aclysmic Variables, ASP Conf. Series, 85, 196, Astron. Soc. 
Pacific., San Francisco 

Wynn, G. A., King, A. R., Horne, K. D., 1997, MNRAS, 286, 436 


APPENDIX A: VALIDITY OF THE LTE 
APPROXIMATION 


LTE requires that the electron, ions and photons are all in 
thermal equilibrium. LTE also requires that detailed balance 
is able to be maintained. Rapid expansion can interfere with 
this. Since recombination rates decrease with density, ion¬ 
ization equilibrium should break down in the low density 
outskirts of the fireball. 

We derive a limiting radius inside which statistical equi¬ 
librium for ionization holds from the condition that the ex¬ 
pansion timescale (£ e xp) be longer than the recombination 
timescale (tree)- Now, 


V _ I Kr 3 r (3 

V 4n r 2 r 3 r 3 H 


(Al) 


and 
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rec ,i 
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OLiTle 


'rii +1 
. rii 
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']} 


-1 


(A2) 


where we have taken the effective recombination timescale 
for the ith state of a given element from Bottorff et al. 
(2000), a; is the recombination coefficient to the ith ion¬ 
ization state and the term in square brackets contains a 
correction factor to the commonly used expression. Here, 
we assume that this term is unity but for detailed checks in 
section Al we evaluate using the calculated densities. 

For an initial estimate, we consider a fireball which is 
composed purely of hydrogen. Equating (Al) and (A2), we 
plot an estimate of th e lim iting ionization equilibrium radius 
for hydrogen in Fig. Al. using parameters M = 10 17 kg, 
ao = 10 s m, To = 18 000 K and a velocity at ao of 
300 km s -1 giving H = 3.0 x 10~ 3 s -1 . We use a radiative 
recombination rate given by 



_ b = 0.7480, T = 3.148 K 

and T2 = 7.036 x 10 5 K ( Verner fe Ferland 1996 ). We can 
estimate the range of times for which the fireball is optically 
thick using a general opacity 
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K, = K,\eT 


\ -3 2 
r V p 


(A4) 


where e is the stimulated emission correction and ki is a con¬ 
stant dependant upon the composition and whether bound- 
free or free-free opacity is being considered. For free-free 
opacity and pure hydrogen composition k\ = 1.34 x 10 52 . 
We parameterize the temperature evolution in the form 
T = To/3 r - Thus, T = 0 for an isothermal model and T = 2 
for adiabatic cooling. Integrating (A4) over a line of sight 
through the centre of the Gaussian density profile we can 
derive an optical depth through the fireball 


Ki tM 2 
Tn 


-r 


Hence the condition for r > 1 becomes 


P'r 


0.35 for T = 0 
0.27 for T = 2 


(A5) 


(A6) 


using the parameters considered above and evaluating at 
3 000 A. We can see that the expanding gas is optically 
thick at early times and that the limiting /3 is relatively 
insensit ive t o the exact parameters given the strong power 
of (3 in ( |A5| ). However, it is most strongly dependant on (at 
least inversely proportional to) the fiducial lengthscale do. 

We can estimate the time for temperature equilibrium 
to be established from the collision frequency between elec- 
trons _and ions using the method of Goldston & Rutherford 

(1995Q 

(AT) 

(AS) 


Tli < CT e iV e > 

2imZ 2 e 4 In A 


127T 2 ejm e 2 (kT e ) 2 

7.3 x 10 -5 SL g" 1 

rp 2 

-L e 

3.3 x 10 8 ^ 3( 5“ 1) e _r ' 2 s -1 . 


(A9) 

(A10) 


where we have evaluated the expression for the case of hy¬ 
drogen using In A = 20 and the same fiducial parameters as 
above. By considering the energy exchanged in a collision 
and integrating over a Maxwellian distribution of electron 
velocities we arrive at a timescale for ions and electrons to 
come into equilibrium 


t e 


3n(2n)2 e%mi(kT e H 

1 

(All) 

n e Z 2 e 4 mi In A 


2.8 x 11 r 6 /3 3(1 ~$ ) e r ' 2 s. 

(A12) 


Electrons collide on a similar frequency but share out there 
energy on a timescale roughly quicker. 

We can estimate the timescale for photons and electrons 
to reach equilibrium in an analogous way using the free-free 
cross-section <Tfj and the opacity outlined above 


n e agc 


p 2 c 

Kie -j 

v 3 T 2 
2.9 x 1(T 


(A13) 

(A14) 


S" 1 (A15) 

where we have evaluated the expression in the uv at 2 000 A. 
This estimate becomes less accurate at low optical depths 
since photons will be able to travel through significant 


Limiting Radii v Time 



Figure Al. Equilibrium limiting radii as a function of time for 
a pure hydrogen fireball. Adiabatic evolutions are shown in solid 
lines, isothermal in dashed. 


changes in density before absorption. This would necessi¬ 
tate an integral over the probability of absorption per unit 
distance for an exact treatment. However, since inward pho¬ 
tons will approach regions of higher density and outward 
lower density the mean density represented by this method 
is sufficiently precise for the use here. Once the expansion 
becomes optically thin, then the assumption of a Planck 
distribution in the derivation of (A15) breaks down and the 
definition of the limiting radius itself becomes more fuzzy as 
photons will interact with electrons throughout the volume. 
When a photon is absorbed it gives up all its energy so that 
we can calculate the timescale to approach equilibrium by 
a similar method to that above. Integrating over frequency 
we arrive at a timescale 


t 


O') 6 


3 c 2 T 2 

167rm.gfvini 

1.8/3 3 ~ie r ' 2 s. 


(A16) 

(A17) 


The limiting radii, where these timescales equal the ex¬ 
pansion timescale are plotted as a function of time in Fig. Al 
alongside the ionization equilibrium limiting radius. The 
lowest of these three boundaries marks the limit for LTE 
to hold. 

Even for the late time behaviour, the minimum photon- 
electron equilibrium radius occurs at R ~ 1.5ao ie. 

2.0cr which contains ~ 95% of the mass. The ionization- 
equilibrium timescale is the most volatile of the plotted 
timescales, since the limiting radius depends on the ion¬ 
ization structure for each element. As such, it should be 
calculated after each simulation using the stored ionization 
history to double check the validity. This radius will also 
evolve differently for each element. 


Al Detailed Check 


Examining the conditions for equality of equations (Al) and 
(A2), we can check the limiting radius for ionization equilib¬ 
rium to hold for each ion included in each of the simulations. 
We plot results, for all the included ions, for the Popula¬ 


tion II adiabatic and isothermal models in Figs. A2 and A3 
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1'onizatian Equilibrium Radius v Time — Adiabatic Population II Mode 


onization Equilibrium Radius v Time - Isothermal Population II Mode 




Figure A2. Ionization equilibrium limit of the simulations for a 
Population II, adiabatically cooling fireball for all ions included 
with non-zero densities. 


Figure A3. Ionization equilibrium limit of the simulations for a 
Population II, isothermal fireball for all ions included with non¬ 
zero densities. 


We use routines to calculate the recombination coefficient 


events bv Verner (Verner & Ferland 1996; Pequij 

?not, Pe- 

tit jean 
Van Stf 
dini & 

& Boisson 1991; Arnaud & Raymond 1992 

Shull & 

enberg 1982; Landini & Monsignori Fossi 1990|; [Lan- 

Monsignori Fossi 1991), dielectronic (Mazzotta et. 


al. 199£), three-body by Cota (1987) as used in Cloudy and 
charge transfer recombination events by Kingdon & Ferland 
(1996). We can see that, in general, there is a slow evolution 
of the limiting radius with time. The quantised nature of the 
grid on which the ionization states, and thus recombination 
timescales, were calculated is apparent (A77 = 0.1). The dra¬ 
matic switching that sometimes occurs results from changes 
in the ionization structure between one spectrum and the 
next. However, since we only require enough spectra to gen¬ 
erate a clear lightcurve, the lack of temporal resolution for 
these changes is not of significant concern. 

Except for late times, both of these evolutions give an 
ionization equilibrium limit rj <; 3.0 ie. 4.2 ct of the density 
distribution. As a result we can be assured that, under our 
assumptions, only a small fraction of the mass of the fire¬ 
ball has non-LTE composition and our LTE treatment is 
reasonable. 


This paper has been produced using the Royal Astronomical 
Society/Blackwell Science ETgX style file. 
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